SUMS OF SQUARES AND ORTHOGONAL INTEGRAL VECTORS 
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Abstract. Two vectors in Z 3 are called twins if they are orthogonal and have 
the same length. The paper describes twin pairs using cubic lattices, and counts 
the number of twin pairs with a given length. Integers M with the property that 
each integral vector with length \J~M has a twin are characterized modulo a famous 
conjecture in number theory. The main tool is the decomposition theory of Hurwitz 
integral quaternions, but elementary, geometric proofs are provided for most of the 
results as well. 



1. Introduction and main results 

An icube in Z n of dimension k is a sequence (i>i, . . . , Vk) of k nonzero vectors in Z n 
that are pairwise orthogonal and have the same length. The subgroup generated 
by Vx,...,Vk is called the corresponding cubic lattice. The common length of the 
vectors Vi is denoted by \\vi\\, and is called the edge length of the icube. By the norm 
of Vi we shall mean ||fj|| 2 (a similar convention is used also for Gaussian integers and 
quaternions). Thus the norm of an integral vector is an integer, but the length is not 
necessarily. 

In this paper we investigate how icubes can be constructed, counted, and extended. 
We mainly consider the case n = 3 and k = 2, but shall make some easy observa- 
tions concerning higher dimensional icubes as well, and formulate the corresponding 
problems in Section 

Before listing the motivating problems and the results in more detail, let us make 
a remark about our methods. Many theorems will be obtained either by elementary 
calculations (Sections El E]) or by geometrical arguments (Section [7]). On the other 
hand, all the results in the paper can be proved using the number theory of Hurwitz 
integral quaternions (see Section [2] for definitions concerning quaternions). These 
arguments occupy Sections HI [5] and [6j We made an effort to keep these sections 
"clean" in the sense that we do not refer here to the proofs in the other sections. 
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Except for some side remarks, we never use the deeper methods of number theory, 
like quadratic extensions, modular forms, analytic tools. However, we believe that 
analytic methods will be useful in dealing with the problems in Section [9j 

Proposition 1.1. Let (v\, . . . ,v n ) be an n- dimensional icube in Z". If n is odd, then 
its edge length is an integer. 

Proof. Let d denote this length. The volume of the cube is d n , which is an integer, 
since it is the determinant of the integer matrix (vi, . . . , v n ). But d 2 is also an integer, 
since the vectors have integer components. Hence d n ~ l is an integer, too. Therefore 
d = d n /d n ~ 1 is rational, and so it is an integer. □ 

Three dimensional icubes in Z 3 have been described by A. Sarkozy |Sar61] . To for- 
mulate his main result, we go back to a construction discovered by Euler. The follow- 
ing well-known facts show how to obtain rotations in M 3 . A rotation is an orientation- 
preserving orthogonal map, that is, an element of the group SO (3). Throughout the 
paper we identify v = (v\, t>2, V3) G M 3 with the pure quaternion V(v) = Vii+V2j+Vzk. 

Theorem 1.2. For a quaternion a define the map R(a) : j3 1— > af3a on the set of 
pure quaternions (3. Then a(3a is also a pure quaternion, so R(a) can be considered 
as a map R 3 — > E 3 . 

(1) If the norm of a is 1 (that is, if a is a unit quaternion), then R(a) is a 
rotation. (In this case R(a) is conjugation by a in the sense of group theory). 

(2) Every rotation is of the form R(a), and this a is unique up to sign. 

(3) Let a = m + ni + pj + qk whose norm is d = \a\ 2 = m? + n 2 + p 2 + q 2 . Then 
R(a) is a dilated rotation, where the dilation coefficient is d. The matrix of 
R(a) in the standard basis (i, j, k) is denoted by E(a). It is the Euler matrix 



Thus the columns of the Euler matrix are pairwise orthogonal, and have common 
length d. We are interested in this matrix when its entries are integers. Call such a 
matrix primitive if the greatest common divisor of its nine entries is 1. Similarly, an 
icube (or a single integral vector) is primitive if the nk entries are relatively prime. 

Theorem 1.3 (Sarkozy, [Sar61] ). Ifm,n,p,q G Z, then E(m + ni+pj + qk) is prim- 
itive if and only if gcd(m, n,p, q) = 1 and d is odd. Every primitive 3- dimensional 
icube in Z 3 can be obtained from such an Euler matrix by permuting columns and 
changing the sign of the third column if necessary. 

This theorem is analyzed in Section [3] and in Corollary 15.121 Sarkozy went on to 
count all icubes with a given edge length d. His proof uses Jacobi's classical result on 
the number of decompositions of integers to the sum of four squares (Theorem 12. 151) . 




—2mq + 2np 
m 2 —n 2 + p 2 — q 2 
2mn + 2pq 
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Now we turn our attention to the question of extension. The above result makes 
it easy to answer the following question: which 1-dimensional icubes (that is, which 
vectors in Z 3 ) can be extended to a 3- dimensional icube? It turns out that the trivial 
necessary condition given by Proposition 11.11 is sufficient. 

Theorem 1.4. A vector in Z 3 is contained in an icube iff its length is an integer. 

Proof. Let u = (a, b, c) be a primitive integral vector whose length d is an integer, so 
a 2 + b 2 + c 2 = d 2 . We may assume that a is odd. It has been known since at least 1915 
(see |Carl5j ) that in this case there exist m, n,p, q G Z such that u is exactly the first 
column of the corresponding Euler matrix. Thus the columns of this matrix extend 
u to the desired icube. 

If a; is a non-primitive vector of integer length, then it can be written uniquely 
as gu, where g G Z and u G Z 3 is primitive. Then the length of u is also an integer, 
so it extends to an icube (u, v,w). Thus (gu, gv, gw) extends x = gu. □ 

Here we referred to the Euler parameterization of Pythagorean quadruples. The- 
orem 18.11 gives the detailed statement with a constructive proof. Theorem 14.21 also 
yields Theorem II .41 using quaternions (see Remark 14.31) . A geometrical proof is given 
in Theorem 17.41 and Corollary 17.51 When u is primitive, the cubic lattice generated 
by any icube containing u is always the same (see Theorem 11.51) . 

The next question is this: which 2-dimensional icubes in Z 3 can be extended to a 
3-dimensional icube? Again the necessary condition that the length be an integer is 
sufficient (Corollary 15.111) . and there is an elementary proof (Proposition 17. ip . 

Having surveyed the investigation of 3-dimensional icubes, we now turn our atten- 
tion to the 2-dimensional case. For brevity, we call a 2-dimensional icube in Z 3 a 
twin pair, and from now on by an icube we shall always mean a 3-dimensional icube 
in Z 3 . The theorems described below are our new results. The essence of them is 
that we understand vectors and twin pairs by putting them into large cubic lattices. 

Theorem 1.5. Let x G Z 3 have norm nm 2 , where n is squarefree. Then there exists 
an icube (u, v, w) with edge length m such that the corresponding cubic sublattice 
contains x. If x is primitive, then this cubic lattice is unique, and is given by an 
Euler matrix E(a) for a quaternion a with integer coefficients. 

The existence part of this result follows from Theorem 14.21 see Remark 14.31 The 
uniqueness part is proved at the end of Section [51 but it is also a consequence of 
Corollary 13.91 and Theorem 14.21 See Corollary 17.51 for an elementary argument in the 
case when x has integer length. 

If (u,v,w) is an icube and a, b G Z, then (av + bw,—bv + aw) is a twin pair. 
Theorem 15.41 says that we get all twin pairs this way. To count all twin pairs, the 
corresponding cubic lattice should be made unique. This is achieved in the same 
theorem by making the cubic lattice as large as possible. However, not necessarily as 
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large as in Theorem 11.51 above. The difficulty is with non-primitive vectors, because 
there is no trivial reduction to the primitive case. For example, 3(8,-10,9) and 
7(4, 5, 2) are twins, and this is explained by the cubic lattice (u, v, w) = E(2i+j +4k) 
with a = 2 and 6=1. Theorem 14.61 tells us how a vector in a cubic lattice can be 
divisible by a prime "unexpectedly" . Theorem 15.41 describes, using the language of 
quaternions, how large this common cubic lattice really is for a given pair of twins. 
As an application, we count all twin pairs with given norm in Theorem 15.101 

The question of extension is more difficult to handle. A consequence of this count- 
ing result is that the common norm of twins is always the sum of two squares (Proposi- 
tion [723 gives an elementary proof). However, the converse is not true, as the example 
of (2, 2, 3) shows: its norm is 17 = l 2 + 4 2 , but it does not have a twin. The case of 
primitive vectors is described as follows (the proof is at the end of Section [5]). 

Corollary 1.6. Using the notation of Theorem \1.5\ suppose that x is primitive and 
x = au + bv + cw. 

(1) If none of a, b, c is zero, then x does not have a twin. 

(2) If exactly one of a, b, c is zero, then x has exactly two twins. If, say, a = 0, 
then these are cv — bw and its negative. 

(3) // two of a, b, c is zero, then x has exactly four twins, and so is contained in 
a unique icube. If, say, a = b = 0, then c = ±1, n = 1, and the twins of x 
are ±u and ±v. This case happens exactly when the norm of x is a square. 

Definition 1.7. A positive integer M is called twin-complete if every vector in Z 3 
with norm M has a twin, and there is such a vector (that is, M is not of the form 
4 n (8A; + 7)). 

Theorem 1.8. A positive integer is twin-complete iff its squarefree part is twin- 
complete. A positive squarefree integer is twin-complete if and only if it can be written 
as a sum of at most two positive squares, but not as a sum of three positive squares. 

We give a complete list of twin-complete numbers modulo the following conjecture. 

Conjecture 1.9. The complete list of those squarefree numbers that can be written 
as a sum of at most two positive squares, but not as a sum of three positive squares 
is the following: {1, 2, 5, 10, 13, 37, 58, 85, 130}. 

Corollary 1.10. The numbers m 2 , 2m 2 , 5m 2 , 10m 2 , 13m 2 , 37m 2 , 58m 2 , 85m 2 , 
130m 2 are twin- complete for every integer m > 0. // Conjecture \1.9l holds, then there 
are no other twin-complete numbers. □ 

The proof of Theorem 11.81 is in Section [6l where many known results concerning 
this conjecture are surveyed. The squarefree numbers in question form a subset of 
Euler's numeri idonei, therefore at most one number can be absent from the list 
above. If such an integer does exist, it must exceed 2 x 10 11 | Wei 73] and if it is even, 
the Generalized Riemann Hypothesis is false [BC00J. 
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2. Integral quaternions 

In this section we list some basic results about the number theory of the ring of G 
of Gaussian integers, and of the ring E of Hurwitz integral quaternions. We assume 
that the reader is familiar with the fact that G has unique factorization, and with 
the description of Gauss primes. We shall need the following observation. 

Lemma 2.1. If a is complex number with rational real and imaginary part and 
a 2 G G, then a G G. 

Proof. Write a = 0/d, where (3 G G and d G Z. Then d 2 | (3 2 among Gaussian 
integers, hence d \ /3, so a = (3/d G G. Here we used unique factorization in G. □ 

Next we review the well-known theorems about decompositions to sums of squares. 

Theorem 2.2 (Fermat). A positive integer d is the sum of two (integer) squares if 
and only if in the prime factorization of d every prime of the form 4k — 1 occurs an 
even number of times. The number of decompositions is 4]T(/ij + 1), where the fij 
are the exponents of the primes of the form 4k + 1 in the canonical form of d. 

This can be proved easily using Gaussian integers. 

Theorem 2.3 (Lagrange). A positive integer d is not the sum of three squares if and 
only if d is of the form 4 n (8k + 7) with some non-negative integers n and k. 

Here it is harder to find the number of decompositions, but the problem can be 
reduced to the case when d is squarefree. The corresponding formula is given by (29) 
in |Pal40] . We include it as Corollary 14.91 because it falls out from the results in 
Section HI 

Theorem 2.4 (Lagrange). Every positive integer is the sum of four squares. 

There is an explicit formula for the number of solutions. We shall state it in 
Theorem 12. 15} because it is related to the number of irreducible quaternions. 

Now we look at the number theory of integral quaternions. The general reference 
is |H19] . The division ring of all quaternions (with real coefficients) is denoted by H. 

Definition 2.5. Let a = a + bi + cj + dk G H be a quaternion (a, b, c, d G R). Then 

a = a — hi — cj — dk 

is the conjugate of a, 

N(a) = a 2 + b 2 + c 2 + d 2 = aa = aa 

is the norm of a, and 

tr(a) = 2a = a + a 

is the trace of a. Therefore 

a 2 - tr(a)a + N(a) =0. 
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A quaternion is pure if its trace is zero. We define the special element 

1+i+j+k 
a= 2 • 

Proposition 2.6 ( [HW79j ). We have a/3 = 0a and N(af3) = N(or) N(/5) for every 
a, (3 el We have N(cr) = f and a 2 = a — 1 . Conjugating by a induces a cyclic 
permutation on {i,j,k} (see Section^ for more details). 

Definition 2.7. Quaternions of the form 

aa + hi + cj + dk (a, b,c,d G Z) 

are called integral quaternions, or Hurwitz integral quaternions; they form a ring E. 
Quaternions of the form 

a + bi + cj + dk (a, b,c,d G Z) 

shall be called quaternions with integer coefficients, or Lipschitz integral quaternions. 
Such a quaternion is called primitive if gcd(a, 6, c, d) = 1. 

Proposition 2.8. A quaternion a = a + bi + cj + dk is Hurwitz if and only if either 
a, b, c, d are all integers, or all of them is the half of an odd integer. If a is such, then 
N(a), tr(a) are integers. A pure integral quaternion has integer coefficients, hence 
the Euler matrix E(a), whose columns are aia, aja, aka, has integer entries. 

We shall be working exclusively in the ring E of Hurwitz integral quaternions, and 
use the divisibility symbol | to denote divisibility on the left in E. The rings Z and G 
are subrings of E, and it is easy to see that here | retains its usual meaning. 

Proposition 2.9. An integral quaternion is a unit (that is, divides every element 
o/E on the left) if and only if its norm is 1. These are exactly the 24 elements 

±l±i±j±k 
±1, ±h ±3, ±k, , 

which form a group under multiplication. Every integral quaternion has a left asso- 
ciate that has integer coefficients ([}IW79\, p. 305). 

Theorem 2.10. The ring E is right Euclidean: for every a, (3 G E with (3^0 there 
exist u, p G E such that a = f3u + p and N(p) < N(/3) (see Theorem 373 of |HW79j ). 

The ring E is also left Euclidean. In fact, a > a is an isomorphism between E and 
its dual, so every assertion that we prove for E holds also if we replace "left" with 
"right" and vice versa. 

As E is left Euclidean, every element can be written as a product of irreducible 
quaternions. This decomposition is unique in a certain sense. For us, it will be 
sufficient to use the following two lemmas for uniqueness. 
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Lemma 2.11. Suppose that a G E and p G Z is a prime such that p | N(a) but 
p does not divide a. Then a can be written as ira' where N(tt) = p, and this tt is 
uniquely determined up to right association. 

An element tt is a left divisor of a with norm p if and only if tt is the generator of 
the right ideal (a,p) r . 

Proof. The ring E is right Euclidean, so the right ideal R = (a,p) r is generated by 
an element tt. Since a,p G R we have that tt divides both a and p on the left. Write 
p = ttt . Then N(7r) N(r) = p 2 . If N(tt) = p 2 , then r is a unit, hence p is an associate 
of ii, and therefore p divides a, contradicting our assumption. If N(tt) = 1, then 
ii is a unit, so ari + pr 2 = 1 for some n, r 2 G E. Thus aari + apr 2 = a showing 
that p | a. Hence again p divides a, which is a contradiction. The only remaining 
possibility is that N(7r) = p (implying that tt is an irreducible element of E). Thus 
7r is a left divisor of a with norm p. 

Conversely, if tt\ is a left divisor of a with norm p, then a and p = ttiWi are in (7Ti) r , 
so (n) r = (a,p) r C (ni) r . Thus tt 1 divides tt on the right, and as they have the same 
norm, they are right associates, implying also that (ni) r = (n) r = (a,p) r . □ 

Lemma 2.12. Suppose that 6, rj, tt G E such that N(tt) = p is a prime in Z. If tt | 9, 
p | 8r] but p does not divide 9, then tt \ 7]. 

Proof. By Lemma 12.111 (p,6) r = (Tr) r , that is, tt = Bt\ + pr 2 for some Ti,t 2 G E. 
Hence Ti] = r{6r] + pf^r], and as p divides 8rj we get that p | Tfrj. Using p = tttt this 
shows that tt \ rj. □ 

Theorem 2.13. An integral quaternion is irreducible in the ring E if and only if its 
norm is a prime in Z (see Theorem 377 of |HW79j ). The only elements o/E whose 
norm is 2 are A = 1 + i and its left associates. If p > 2 is a prime in Z, then there 
exist exactly 24(p + 1) integral quaternions whose norm is p. 

Corollary 2.14. The number of integral quaternions with norm n is 24 times the 
sum of positive, odd divisors of n. 

Theorem 2.15 (Jacobi). Let a(s) denote the sum of positive divisors of any integer s, 
and for every positive integer n, let N be the number of solutions of the Diophantine 
equation 

x 2 + y 2 + u 2 + v 2 = n (x, y, u, v G Z) . 

Ifn is odd, then N = 8a(n). If n = 2 t m, where m is odd andt > 0, then N = 24a(m) 
(see Theorem 386 of |HW79j ). 

Lemma 2.16. Let p G Z be a prime and t > 0. Suppose that tt\ G E is fixed, and 
has norm p. Consider all integer quaternions a such that N(a) = p and air\ is not 
divisible by p. Then the number of such a is 24p £ . 
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Proof. We do induction on I. If I = 0, then the statement is trivial, since the number 
of units is 24. Suppose that p does not divide a. The dual of Lemma \2 . 1 1 1 shows that 
a can be written as a 2 7r 2 , with N(7r 2 ) = p, where a 2 is unique up to right association, 
and 7r 2 is unique up to left association. Apply Lemma 12.121 for 9 i— > a, rj i— > tti and 
7r i — ► 7T2- We get that if a7Ti is divisible by p, then 7r 2 and 7fT are left associates. 
Conversely, if 7r 2 and 7fT are left associates, then clearly p | mri. 

By Theorem 12. 13[ the number of elements of norm p up to left association is p+ 1. 
So 7r 2 can be chosen p ways, and by the induction assumption, a 2 can be chosen 
24p e ~ 1 ways for every given 7r 2 . Thus a can be chosen 2Ap e ~ 1 p ways. □ 

3. Integral Euler matrices 

In Theorem 11.31 describing primitive icubes Sarkozy uses Lipschitz integral quater- 
nions. It will be more convenient in this paper to generate cubic lattices by Hurwitz 
quaternions. Our goal in this section is to characterize (in Theorem 13.31) all Euler 
matrices E(a) with integer entries (called integral Euler matrices) in terms of the 
corresponding quaternion a. As Corollary 13.91 Sarkozy 's theorem is obtained. The 
proofs are elementary in the sense that decomposition to irreducible quaternions is 
not used. 

We now explain how to permute the columns of an Euler matrix by changing its 
generating quaternion. By Theorem 11.21 E(a) is the matrix of R(a) : f3 \— > a/3 a, 
hence E(ae) = E(a)E(e). The map R(a) is always orientation-preserving, but the 
map corresponding to an icube (as a matrix) is not necessarily. This problem is 
averted by taking the negative of an odd number of columns. 

Proposition 3.1. Let e be 

(A) a or <7 -1 , where a — (1 +i + j + k)/2. Then R(e) is the rotation o/R 3 about 
the vector i + j + k through an angle o/±120° (thus cyclically permuting the 
three coordinate axes). Therefore E(ae) is obtained from E(a) by applying a 
cyclic permutation to the columns. 

(B) (l±z)/v^2- Then R(e) is the rotation about the unit vector i e Z 3 through an 
angle o/±90° (hence interchanging the other two coordinate axes). Therefore 
E(ae) is obtained from E(a) by switching the last two columns and taking the 
negative of one. A similar statement holds for (1 ± j)/v^2 and (1 ± k)/ ' \[2. 

(C) ±i. Then R(s) is the half turn (that is, 180° rotation) about the unit vector 
i G Z 3 (fixing all coordinate axes). Therefore E{ae) is obtained from E(a) by 
taking the negative of the last two columns. A similar statement holds for ±j 
and ±fc. This transformation is the square of the one described in (B). 

Every non-identical permutation of the columns of E(a) can be obtained by one of 
the above modifications of a, but in case of an odd permutation one of the columns 
changes its sign. One can also change the sign of any two columns. □ 
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Before going further let us record the action of these isometries on IR 3 . 

Proposition 3.2. Denote by H the group of units of E (see Proposition \2.9\) . set 
Q = {±1, ±i, ±j, ±k} and let G be the subgroup of the multiplicative group of M 
generated by H and (1 + i)/y2. Then G contains all the isometries investigated in 
Proposition \3.1[ The group G has 48 elements. 

The element a has order 6. The rotation R(o~) maps rj = ai + bj + ck to ci + aj + bk, 
so it permutes the components cyclically. 

The element (1 + has order 8. The corresponding rotation + i)/V2j 

maps 7] to ai — cj + bk. The square of this rotation is R(i) mapping rj to ai — bj — ck. 

In general, G acts on the set of pure quaternions via the rotations R(p) with p G G. 
The orbit of i] under Q consists oft], —ai — bj + ck, —ai + bj — ck and +ai — bj — ck. 
If we disregard the signs, then every other element of H induces a fixed point free 
permutation on the components ofrj. □ 

The proof is left to the reader. Note that these groups H and G are the so-called 
binary tetrahedral group and binary octahedral group, respectively. 

Theorem 3.3. E(a) is a primitive integral Euler matrix if and only if the non-zero 
quaternion a belongs to one of the following three types. 

(1) a is a primitive Lipschitz integral quaternion with an odd norm. 

(2) a = /3/y/2, where (3 is a primitive Lipschitz integral quaternion such that 
N(/3) = 2 (4), or equivalently: exactly two components of (3 are odd. 

(3) a = (3/2, where (3 is a primitive Lipschitz integral quaternion such that 
N(/3) = 4 (8), or equivalently: all four components of (3 are odd (so a is a 
Hurwitz integral quaternion). 

In all cases N(a) is an odd integer. Each column and each row of E(a) contains 
exactly one odd entry. The number of odd entries in the main diagonal in types (1), 
(2), (3) are 3, 1, 0, respectively. 

The proof will utilize several lemmas. The first one explains the equivalences stated 
within the theorem. 

Lemma 3.4. If (3 is a primitive Lipschitz integral quaternion, then N(/3) cannot be 
divisible by 8. It is congruent to 2 modulo 4 iff (3 has exactly two odd components, 
and is congruent to 4 modulo 8 iff all components of (3 are odd. 

Proof. As (3 is primitive, it has an odd component. Thus the statement follows from 
the fact that the square of an odd integer is congruent to 1 modulo 8. □ 

Lemma 3.5. If E{a) is a primitive integral Euler matrix, then N(a) is an odd integer. 
Each column and each row contains exactly one odd entry. 
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Proof. By Theorem 11.21 the length of every column and row of E(pt) is N(a), which 
is an integer by Proposition 11.11 Suppose that it is an even number. Let (a, b, c) be 
a row of E(a). Then taking a 2 + b 2 + c 2 = N(a) 2 modulo 4 we see that a, b, c are 
even. Thus each entry of E(a) is even, violating the condition that E(a) is primitive. 
Therefore a 2 + b 2 + c 2 = 1 (4) showing that exactly one of a, b, c is odd. □ 

The next lemma contains an elementary proof of the last part of Proposition 12.91 
for right associates. 

Lemma 3.6. If the quaternion a = (m + ni + pj + qk)/2 belongs to class (3) of 
Theorem \3.3\ then either aa or aa' 1 is a quaternion of class (1). 

Proof. Elementary calculation shows that 

m — n — p — q + (m + n + p — q)i 



aa 



aa 



4 

(m — n + p + q)j + (m + n — p + q)k 
+ 4 
m + n + p + q+ (— m + n — p + q)i 



+ 



4 

-m + n + p — q)j + (— m — n + p + q)k 



4 

By Lemma 13.41 all integers m, n, p, q are odd. Therefore, either m — n — p — q oi 
m + n + p + q is divisible by 4. However, all four integer coefficients in the numerator 
of aa are congruent to each other modulo 4, and the same is true for aa -1 , hence one 
of these two quaternions is a Lipschitz integral quaternion with an odd norm N(a). 
The components of this quaternion are obviously relatively prime. □ 

Lemma 3.7. Every quaternion a = (m + ni + pj + qk)/\/2 of type (2) can be 
multiplied on the right by a suitable unit (l + w)/v^2 to transform it to type (1), where 
u e {i,j,k}. 

Proof. By Lemma I3~4l exactly two of m, n, p, q are odd. We may assume that m and n 
have the same parity (dealing with the other two cases is analogous). Thus 

a = + \/2(a + bi + cj + dk) or a = j -j= + \pl{a + bi + cj + dk) 



with a,b,c,d E Z. Obvious calculation shows that the quaternion a{l + i)/y2 = (3 
has integral coefficients and N(/3) = N(a) is odd. □ 

Proof of Theorem \3.3[ Put a = m + ni + pj + kq with real numbers m, n,p, g, and 
assume that the Euler matrix E(a) given in Theorem 11.21 has integral entries and is 
primitive. By Lemma 13.51 N(a) = m 2 + n 2 + p 2 + q 2 is an integer, and the diagonal 
elements of E(a) are integers, too. Taking linear combinations of these quadratic 



SUMS OF SQUARES AND ORTHOGONAL INTEGRAL VECTORS 



11 



forms we get that 4m 2 , 4n 2 , 4p 2 , and 4q 2 are all integers. By adding and subtracting 
symmetric off-diagonal elements we also obtain that all the numbers 4mn, 4mp, Amq, 
4np, 4nq, and 4pq are integers as well. Therefore, the squarefree parts of the non-zero 
numbers among 4m 2 ,4n 2 ,4p 2 and 4q 2 are the same. Denote this common squarefree 
part by r. The quaternion a = m + ni + pj + qk can be written uniquely in the form 

(3.8) a = ^^(a + bi + cj + dk), 

where k G N, a,b,c,d G Z, (a,b,c,d) = 1. Since the matrix E(a) is primitive, 
neither k, nor the squarefree r can have any odd prime divisor, and k (as a power 
of 2) cannot be greater than 2. Hence both k and r are elements of the set {1,2}, 
but k = r = 2 violates primitivity of E(a). Thus, we are left with the cases 

(a) k = 2, r = 1; 
(6) k=l,r = 2; 
(c) k = r = 1. 

These correspond exactly to the cases listed as (1), (2), and (3) in Theorem 13.31 
Lemma 13.51 shows that N(ct) is an odd integer, therefore Lemma 13.41 finishes the 
proof of one implication of the theorem. 

Assume now that the quaternion a is one of the types (1) — (3) in Theorem 13.31 
We want to show that E(a) is a primitive integer matrix. By Lemmas 13.61 and 13.71 
there exists a suitable quaternion e G H with N(e) = 1 such that ae is of type (1), 
and by Proposition 13.11 we see that E(ae) = E(a)E(e) is a primitive integral matrix 
if and only if E(a) is. 

We show that E(ae) is primitive. Since N(ae) is odd by assumption, the entries in 
the main diagonal of E(ae) are odd. By way of contradiction, suppose that an odd 
rational prime t divides all entries of E(ae). Let ae = m + ni + pj + qk. The simple 
calculation preceding (13.81) shows that t divides the numbers 4m 2 ,4n 2 ,4p 2 , and 4g 2 , 
violating primitivity of ae. Thus, E(ae) = E(a)E(e) is indeed a primitive integer 
matrix. 

Finally we show the last statement of the theorem. If a belongs to type (1), then, 
as we saw above, the entries in the main diagonal are odd, while the other entries are 
clearly even. Type (2) quaternions are handled by Lemma 13.71 and these correspond 
to the switch of two columns by (B) of Proposition 13.11 The type (3) case is handled 
by Lemma I3.6[ and yields a cyclic, fixed point free permutation of the columns by 
(A) of Proposition 13.11 This completes the proof of Theorem 13.31 □ 

Corollary 3.9. Consider a primitive icube as the columns of a matrix M. Then 
there exists a Lipschitz integral quaternion a such that by permuting the columns 
of E(a), and changing the sign of the last column if necessary, we get M. 

Proof. Change the sign of the last column iff M is orientation-reversing. The new M 
can be written as M = E(a) for some quaternion a (with real coefficients), by 
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Theorem 11.21 This a has one of the three types of Theorem 13.31 Modify a using 
Proposition 13.11 so that the odd entries move to the main diagonal. Then we get an 
integral Lipschitz quaternion by the last statement of Theorem 13.31 If this transfor- 
mation changes the sign of a column other than what has originally been the third, 
then use Proposition 13.11 to change the sign of two columns back. □ 

4. A REPRESENTATION OF PURE INTEGRAL QUATERNIONS 

In this section we decompose single pure quaternions. Geometrically, this means 
that we find a large cubic lattice that contains the corresponding vector. Alge- 
braically, a cubic lattice is the subgroup of all quaternions 8 = a/3 a, where a is a 
fixed Hurwitz integral quaternion, and j3 runs over all pure integral quaternions (see 
Corollary 15.111) . The generating icube is given by ma, aja, aka. 

The primitive case is easier, and is handled by Theorem 14.21 This already implies 
Theorem 11.41 and the existence statement of Theorem 11.51 Theorem 14.61 explains 
how a vector in a cubic lattice can be divisible by a prime "accidentally" . This will 
be used in the characterization of twin-complete numbers, and is also sufficient to 
obtain a classical result about counting all vectors of a given length (Theorem 14.81 
Corollary S3]). 

The results in this section are closely related to those in |Pal40] . but that paper 
deals with Lipschitz quaternions. 

Lemma 4.1. Let 8 G E be a pure quaternion and p G Z a prime such that p 2 | N(<5) 
but p does not divide 8. Then there exists an element it G E whose norm is p such 
that 8 = it5iT for some 8\ G E. 

Proof. By Lemma [2.111 we get that 8 = 7r5 2 for some tt, 82 G E such that N(tt) = p. 
Then N(5) = pN(5 2 ), hence p divides N(<5 2 ) but p clearly does not divide 8 2 - By the 
dual of Lemma [2.111 we obtain an element ir\ G E with norm p such that 82 = 83^1. 
Hence 5 = ^837^1. Taking conjugates we get 8 = W18371. But 8 is a pure quaternion, 
hence 8 = —8. Therefore 8 is divisible by it and by Wi on the left. By the uniqueness 
statement of Lemma [2. Ill we get that n and 7T1 are right associates. Thus 8 = 1183711 
can be written as ttSiT indeed. □ 

Theorem 4.2. Let 8 G E be a pure quaternion whose norm N(£) = nm 2 . Suppose 
that no integer prime divisor of m divides 8. Then 8 can be written as a(3a for some 
a, (3 G E such that N(a) = m and N(/3) = n. Here a is uniquely determined, that 
is, any two such elements a are right associates of each other, and the corresponding 
elements (3 are group-conjugates of each other via a unit of E. If n = 1, then (3 can 
be chosen freely to be any element of {±i, ±j, ±k}. 

Proof. The existence of a and f3 is easily proved by induction on m: apply Lemma I4TT1 
successively for the prime divisors of m. 
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For the uniqueness assume that 5 = ai(3{ai = a^/^a^. We use induction on m 
again. If m = 1, then aii and a 2 are units, so they are right associates, and the 
unit e = oi^ oil satisfies that e(3ie~ l = 02- If Jn > 1, then let p G Z be a prime 
divisor of m. Apply Lemma [2.111 to get m, 7r2 G E with m | ot\ and H2 \ ol%. Then m 
and 7T2 divide 5 on the left, hence the uniqueness statement of Lemma [2.111 implies 
that 7Ti and 7r 2 are right associates. So if 7r 2 = me, a± = nia^ and «2 = 7r 2 a;4, 
then 5' = a 3 /?i = (ea^fasal. By the induction hypothesis, a 3 and ea 4 are right 
associates. Thus so are a\ = nia 3 and a 2 = mtoti- 

If n — 1, then /3 is a unit in E. Since /3 is a pure quaternion, it is contained in 
{±i,±j,±k}. These six elements are group-conjugates of each other via a unit by 
Proposition 13.21 and therefore by taking a right associate of a we may choose any of 
them to be (3. □ 

Remark 4.3. Theorem 14.21 implies Theorem II A\ and the existence statement of The- 
orem [T3J (The uniqueness part of Theorem II .51 clearly follows from Theorem 14.21 and 
Corollary I3.9| but we give a "pure number-theoretic" proof in Section 0) 

Proof. Let u be a primitive vector and denote by 5 the corresponding pure quaternion. 
Decompose 5 using Theorem 14 .21 as 5 = a/3 a with N(a) = m. Then the cubic lattice 
corresponding to a has edge length m and contains u. This yields the existence 
statement of Theorem 11.51 

If the length if u is an integer, then n = 1 and we may assume that f3 = i. Then 
aj a and aha extend u to an icube, proving Theorem 11.41 in the primitive case. The 
general case obviously follows from this. □ 

Lemma 4.4. Let (3,/3i & E be pure quaternions whose norm is the same number n. 
IfP + fr^O, then we have (/? + /3i)/3(/3 + A)" 1 = fa- 

Proof. The proof is an easy calculation using (3 2 = f3\ = —n. Instead of presenting it, 
we explain this formula geometrically. Since 7 = f3 + (3\ is a nonzero pure quaternion, 
conjugation by 7 acts on M 3 as half turn about the line through 7, which clearly takes 
(3 to /3 X . □ 

Lemma 4.5. Let n,j3 G E such that (3 is a pure quaternion and p = N(tt) > 2 is a 
prime in Z. Then n/Sn^ 1 elf and only if there exists an integer h G Z such that 
7f I h + p. 

Proof. First suppose that W \ h + (3, that is, Wt = h + (3 for some r G E. Then 
prW = 7T7fr7f = nhW + tt(3W = ph + 7T/3 7f . Hence p | ir(3W, which implies that 
7T/37T -1 = (tt/3W)/p is indeed an integral quaternion. 

To prove the converse set (3\ = n(3^~ x and r = (3 + (3\. We can assume that n does 
not divide r on the left, as we now show in this paragraph. Let 

(3 2 = (ztO/^tt)- 1 = iftr 1 , 
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this is still an integral quaternion. It is clearly sufficient to prove that m \ h + P, so 
we can work with in instead of 7r in the argument below. However, if both 7r and in 
are "bad", that is, n \ r = (3 + Pi and also in \ P + f3 2 , then it \ i~ 1 j3i + (3\, and so 
Ti | P — i~ l pi. Put P = ai + bj + ck, then P — i~ x Pi = 2{bj + ck). If jn and kn are 
"bad", too, then n divides 2(ai + ck) and 2{ai + bj) as well. Taking norms we get, 
using N(7r) = p > 2, that p divides a 2 + b 2 , a 2 + c 2 and b 2 + c 2 , so p divides a, b, c, 
and therefore p divides p. Therefore 7f | h + P for h — 0, and so we are done in this 
case. Thus we can indeed assume that n does not divide r on the left. 

Lemma S3 implies that r/?r _1 = Pi = nPn~ l , so r _1 7r centralizes p. Let d = N(r), 
then dr~ lr K = ttx centralizes P as well. The centralizer of P consists of the elements 
r + sP, where r,s G R. This set is closed under conjugation, since P is a pure 
quaternion, and therefore contains fix = tit. Write W t = u + vP, where u and t> are 
real, and P = dp', where d G Z and /3' is primitive. Then 2u and 2v<i are integers. 
(We need the factor 2, because the integral quaternion Wt need not have integer 
coefficients). 

Next we show that p does not divide 2vd. Suppose it does. Taking norms we get 
that N(vf) = p | 4 N(u + vp) = (2m) 2 + (2vd) 2 N (/?'), so p | 2u. Thus either u + vp has 
integer coefficients, which are divisible by p, or 2u + 2t> /3 has odd integer coefficients 
which are divisible by p. Since p ^ 2, we have u' + t>'/3 = (u + v0)jp G E in either 
case. Then Wt = u + vP = p{u' + v'P) = Wir(u r +v'P), so r = ir(u'+v r P), contradicting 
our assumption that it does not divide r on the left. Therefore we indeed have that 
p does not divide 2vd. 

Let x, y be integers such that (2vd)x + yp = 1. Then 

7f I 2x(u + vdp') = x(2u) + (1 - yp)p' . 
Since W \ p we get that 7f | x(2u) + /?' so the integer h = x(2u)d works. □ 

Theorem 4.6. Suppose that a,P G E siic/i t/iat P is a pure quaternion and p G Z is 
a prime. Then p | aP a if and only if one of the following cases hold. 

(1) p divides a or p. 

(2) p = 2, which does not divide a,P but divides N(a). 

(3) p > 2, which does not divide a,P but divides N(a) and there exists a right 
divisor it of a with norm p and an integer h G Z such that W | h + p. 

In particular, every prime divisor of aPa divides either J3 or N(a). 

Proof. If (1) holds, then clearly p | aP a. If (2) holds, then the dual of Lemma [2.111 
shows that a is right divisible by 1 + i (since this is the only element in E of norm 2 
up to left association), and Proposition 13.21 yields that (1 +i)P 1 + i is divisible by 2. 
Finally if (3) holds then p | 7r/?7f by Lemma 14.51 This proves one direction of the 
theorem. 

Now assume that p \ a Pa but a and P are not divisible by p. If p does not divide 
N(a), then we have p \ a(aPa)a = N(a) 2 /3, hence p \ P, a contradiction. Hence we 
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may assume that we are in case (3), that is, p > 2 and p | N(a). We proceed by 
induction on N(a). By the dual of Lemma [2.111 a = a\K for some ai,7T G E with 
N(7r) = p. We show that Pi = ir(3w is divisible by p. Then we are clearly done by 
Lemma 14.51 

Suppose to get a contradiction that Pi is not divisible by p. Apply the induction 
hypothesis to ai/3ia7 (which is equal to a/3 a). We must be in case (3), since p > 2 
and p does not divide a,\ and P\. Therefore there exists a right divisor tti of ot\ 
of norm p and an integer h\ such that 7fT | hi + Pi. Taking norms we see that 
p = N(7fT) | N(h +P) = hl + N(Pi). But N(p x ) = p 2 N(/3) is divisible by p, so 
p \ hi, and therefore W[ \ Pi. Since /3i is not divisible by p, the uniqueness part of 
Lemma 12.111 shows that Wi and it are right associates. But then a is divisible on the 
right by 7Ti7fi = p, contradicting our assumptions. □ 

Proposition 4.7. Suppose that (3 G E is a pure quaternion and p is a prime not 
dividing p. Denote by e the number of different quaternions of the form ePe~ l , where 
e runs over the units o/E. Consider all quaternions a whose norm is p £ with some 
fixed £ > 0. If p > 2, then the number of quaternions of the form aft a that are not 
divisible by p is 

(1) ep e ifp\N(P); 

(2) e(p e — p^ 1 ) if — N(P) is a quadratic residue modp; 

(3) e(p £ + p^~ 1 ) otherwise. 

If p = 2, then this number is 0. 

Proof. If p = 2 (and i > 0), then Theorem 14.61 shows that aP~a is divisible by p, so 
suppose that p is odd. Call a pair (ck 1; Pi) "good" , if ai G E with N(ai) = p e , there is 
a unit e G E such that /?i = e^fie, and p does not divide cuiPi aT- By the uniqueness 
part of Theorem 14.21 every element aiPi ai is given by exactly 24 pairs. Therefore it 
is sufficient to count the good pairs for any given p x . 

Let (ai, Pi) be a good pair, then clearly ati is not divisible by p, so we can write 
a.i = a2 7r 2 by Lemma I2.11[ where tt 2 of norm p is uniquely determined up to left 
association. Theorem 14.61 shows that ceiPiOi is divisible by p if and only if TfJ divides 
h + Pi for some integer h (assuming that ai is not divisible by p). Let us count the 
number of such quaternions 7r 2 . 

Clearly, W5\h + p 1 implies N(7r 2 ) = p | N(/i + Pi) = h 2 + N(/3i). This means that 
either p | N(/?i) or — N(Pi) is a quadratic residue mod p. Hence in case (3) above 
there is no such 7r 2 . Since p does not divide Pi, it does not divide h + Pi. Thus by 
the uniqueness statement of Lemma 12.111 there is exactly one left divisor up to 
right association with norm p of any given h + Pi for which p | h 2 + N(Pi). Thus 7r 2 
is unique up to left association. Clearly, the numbers hi and hi yield the same if 
and only if hi = /i 2 (p). If p | N(/?i), then fa = is the only possibility, this yields 
one "bad" value for 7r 2 up to left association. Otherwise there are exactly two values 
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1 < h < p — 1 such that p | h 2 + N(/?i) (since p is an odd prime, assuming of course 
that — N(/?i) is a quadratic residue modp). So in this case there are two "bad" values 
for 7r 2 up to left association. 

By Theorem 12 . 1 3 1 1 he number of possible choices for n 2 is p+1 up to left association. 
Thus for every given (5\ the number of good values for tc 2 is p, p — 1 and p + 1, 
respectively, corresponding to cases (1), (2) and (3) in the claim. If tt 2 is fixed, then 
the number of choices for a 2 , so that ot\ is not divisible by p, is 24j/ _1 by Lemma f2. 161 
Since the number of possible fi\ is e, we get the result. □ 

In the theorem below (— n/p) denotes the Legendre symbol (which is defined to 
be if p | n). 

Theorem 4.8. Suppose that m,n > 1 are integers and n is squarefree. If m is odd, 
then the number p(nm 2 ) of primitive vectors (x,y,z) whose norm is nm 2 is 



where p runs over the prime powers in the canonical form of m. If m is even, then 
p{nm 2 ) = 0. 

Proof. We do induction on the number of prime divisors of m. Let m = j/mj, where 
p does not divide mj. Theorem 14. 21 implies that the pure quaternion 5 = xi + yj + zk 
corresponding to (x,y, z) can be represented as a/3 57, where N(a) = p l . 

Theorem 14.61 shows that if (3 is primitive, then the only possible prime divisor of 
a(3 a is p, and if p = 2 and £ > 0, then a(3a is not primitive, because it is divisible 
by 2. Thus if p > 2, then a/3 a is primitive if and only if j3 is primitive and a/3 a is not 
divisible by p. Then the formula for p(nm 2 ) follows clearly from Proposition 14. 71 □ 

The following is a well-known formula (see (29) in |Pal40] ) . and it follows with 
some effort from Theorem 14.81 

Corollary 4.9. Suppose that m, n > 1 are integers and n is squarefree. The number 
s(nm 2 ) of all vectors of norm nm 2 is 



where p runs over the odd prime powers in the canonical form of m and a(s) denotes 
the sum of positive divisors of any integer s. 



Theorem 15.41 is our main characterization of twins. In Theorem 15.101 we count 
twin pairs with a given norm. Finally we deal with the problem of extension. In 
Corollary 15.111 we show that each pair of twins whose length is an integer extends to 
an icube. Then, at the end of the section, we prove Theorem 11.51 and Corollary 11.61 




s{nm 2 ) = s(n) Y\ (°"G° £ ) — (~ n /p) cr (p e X ) 
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First we translate everything to the language of quaternions. Recall that every 
vector v = (vi, v 2 , v 3 ) G IR 3 is identified by the pure quaternion V(v) = Vii + v 2 j + v 3 k. 

Proposition 5.1. Two vectors v and w in 1? are twins if and only if 9 = V(v) and 
rj = V(w) satisfy the following conditions. 

(1) 9 and rj are nonzero pure quaternions; 

(2) N(0) = Nfa); 

(3) 9r) is also a pure quaternion; 

(4) 9 and rj have integer coefficients. 

We call a pair of such quaternions (9, n) a twin pair. 

Proof. It is easy to verify that the real part of V(v)V(w) equals the negative of the 
dot product of v and w, and the pure quaternion part of V(v)V(w) corresponds to 
the cross product of v and w. □ 

Next we translate the construction of twin pairs given in the Introduction. Denote 
by (u, v, w) the columns of an Euler matrix given by a G E (the ring of Hurwitz inte- 
gral quaternions, see Definition 12. 71) . By Proposition 12.81 N(q) and the components of 
(u, v, w) are integers. Let z = a + hi G G (the ring of Gaussian integers). Then E(a) 
maps i to V(u), j to V(v), k to V{w). It maps the twin quaternions zj = aj + bk 
and zk = ak — bj to azj a and azk 57, respectively, which therefore correspond to the 
twin pair (av + bw, —bv + aw). 

Definition 5.2. We say that (9, if) is parameterized by the pair (a,z) G H x C if 
9 = azj a and i] = azk a. Two pairs in (ai, z\) G HI x C and (a 2 , z 2 ) G EI x C are 
equivalent if they parameterize the same (9,7]), that is, if a\Zxfa{ = a 2 z 2 ja 2 and 
a\Z\ka~i = a 2 z 2 ka 2 . 

Proposition 5.3. Suppose that (9, rj) is parameterized by a pair (a, z) G EI x C. Then 
9rj = N(a) N(z)ai a, so 9 and r\ satisfy (1) — (3) of Proposition [3T71 If in addition 
a G E and zeG, then 9 and rj have integer coefficients, hence they are twins. 

Proof. Note that a~ l = a/ N(a) and x i— > axa~ x is an automorphism of the di- 
vision ring EI. Therefore 9rj = N(a)azjzka. But zjzk = zjzj~ l jk = zzi, so 
9rj = N(a) N(z)aia, which is pure. □ 

Theorem 5.4. The characterization of twin quaternions is the following. 

(1) The quaternions 9 and rj are twins if and only if (9,rj) is parameterized by a 
pair inExG (whose components are nonzero). 

(2) Every pair inExG is equivalent to a pair where the second component is 
squarefree in G. 

(3) Let (a±, z{), (a 2 , z 2 ) G E x G be such that both Z\ and z 2 are squarefree. Then 
these pairs are equivalent if and only if there exists a unit p G G (that is, an 
element o/{±l,±i}J such that a 2 = a%p and z 2 = p 2 zi. 
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(4) The length of the twins 9 and rj is an integer if and only if in the parameter- 
ization (a,z) of (9,rj) where z is squarefree, the number z G G is either real 
or pure imaginary. 

The condition that z is squarefree expresses the fact that the cubic lattice given 
by a is as large as possible. We prove this theorem through a series of assertions. 

Lemma 5.5. Suppose that 9, r) G E such that 9, i] and 9r\ are pure quaternions. Let 
p G Z be a prime such that p divides N(#) and N(r/). Then p | 9r\. 

Proof. Suppose that p does not divide 9rj. By Lemma [2.111 and Lemma [4.11 we have 
9 = 71"!$!, r] = "i]xK2 and 9r] = tt5iW, where tt, tt\, TT2 have norm p. By the uniqueness 
part of Lemma 12.111 tt, 7i"i and W2 are right associates. However, 9, rj and 9rj are 
pure quaternions, so 9i] = —9i] = —Tj9 = —i]9 = —r\Yi\2^\9\ is divisible by p, a 
contradiction. □ 

Lemma 5.6. Suppose that 9, 77 G E such that 9, 77 and 9r] are pure quaternions. Let 
p G Z be a prime such that p 2 divides N(#) and N(?/) ; but p does not divide 9 and rj. 
Then there exist elements tt, 9i, 771 G E such that N(7r) = p, 9 = tt9{t\ and i] = nrjiW. 

Proof. By Lemma 15.51 p \ S = 9rj. Using Lemma 14.11 we can write 9 = tt9{ti and 
rj = 7Tii]iWi, where N(7r) = N(7Ti) = p. Since 9 is pure, 5 = 9rj = —9i], so Lemma \2. 121 
shows that tt \ rj. Applying the uniqueness part of Lemma 12.111 to 77 we obtain that 
tt and iTi are right associates. □ 

Lemma 5.7. Suppose that 9, rj G E such that 9, rj and 9r] are pure quaternions. Let 
p G Z be a prime such that p 2 divides N(#) but p does not divide 9. Then there exist 
elements tt, 9±, r\\ G E such that N(7r) = p, 9 = t\9{t\ and pi] = t\i]{tx . 

Proof. Again write 9 = tt9iW. Suppose first that tt | 77, that is, 77 = ttt]2 for some 
i]2 G E. Then pi] = tt(i]2Tt)T, so we are done in this case. So we can assume that 
tt does not divide 77. Since 9 is pure, 5 = 9rj = —9i], so Lemma 12.121 shows that 
p does not divide 5. By Lemma 14.11 and the uniqueness part of Lemma I2.11[ we 
have 5 = ttSiW. Then Tr9iWri = 5 = Tx5\Tf implies that 9iWi]tt = 5itttt = p8\. Hence 
p I 7f777r 9\. Lemma [2.121 shows that either Tf | 9\ or p \ Wi]n. The first case is impossible 
because then 9 would be divisible by p. Let Wi]tt = pi]i, then pi] = icqin. □ 

Proposition 5.8. Every pair (9, if) of twins is parameterized by a pair (a, z) G E x G. 

Proof. Let M = N(#) = N(r7). Consider all representations of the form 9 = a9{a and 
77 = ai]ia, where a, #1,7/1 G E. Clearly, 9\ and 7/1 are twins as well. There exists such 
a representation (with a = 1), so there is one where N(a) is as large as possible. We 
present reductions to increase N(a). 

If p is a prime such that p 2 divides both N(#i) and N(i]i), but p divides neither 9\ 
nor 7/1, then Lemma [5.61 allows us to replace a with an. If this p does not divide 9\ 



SUMS OF SQUARES AND ORTHOGONAL INTEGRAL VECTORS 



19 



but divides rji, then write 771 = p e r]' such that 7/ is not divisible by p, and apply 
Lemma 15.71 to 9\ and rj'. Then again we obtain a suitable tc (by using up a factor 
of p out of p e ). 

If none of these reductions can be performed anymore, then 9\ = d6 2 and 771 = dt] 2 
for some d G Z such that N(# 2 ) = N(r/2) is squarefree. By Lemma [5.51 every integer 
prime divisor of N^) divides 9 2 7] 2 . Hence the squarefree number N(r/2) divides ^2^2, 
whose norm is N^) 2 . Therefore Q 2 i] 2 = N(r/2)£, where e is a unit of E. We may 
assume that e = i by the argument in the last paragraph of the proof of Theorem 14.21 
Thus 772 = 6 2 i , and as #2 and 772 are pure quaternions, 6 2 = Z\j for some z 1 G G. 
Then = a(dzi)ja and 77 = a(c£zi)A;a, proving the claim. □ 

Lemma 5.9. Let (a,z) G H x C and suppose that z = s 2 t for some s,t G C. Tnen 
tae pairs (a,z) and (as, i) are equivalent. 

Proof. Since jsj^ 1 = s for every s G C, we have that = stjs and similarly 
zA; = stks. Therefore azj a = (as)tj(as) and azka = (as)tk(as). □ 

This lemma immediately implies (2) of Theorem 15.41 and (1) has been proved in 
Proposition 15.81 We now proceed to prove (3). If a suitable p in (3) exists, then 
(ai, zi) is equivalent to (a±p, z%) = («2, z 2 ) by Lemma [5791 So we only have to prove 
that if (ai, zi) and («2, ^2) are equivalent pairs in E x G, then a suitable p exits. 

Choose elements s r 6C satisfying si = z r . Lemma [5791 shows that (a r , z r ) is equiv- 
alent to (a r s r , 1). From (aiSi)j(aTsi) = («2S2)j(«2S2) we get that (a2S2)~ 1 (aiSi) 
centralizes j, and it centralizes k, too, by the same calculation. Since the ele- 
ments of H centralizing both j and k are exactly the real numbers, we get that 
t = («2S2) _1 (aiSi) is contained in K. This can be written as taf 1 *^ = S1S2 . From 
(aiSi)j(oiSi) = (a 2 s 2 )j(a2S2) we get that N(a!) 2 N(si) 2 = N(a 2 ) 2 N(s 2 ) 2 , hence 
N(t) = 1, and thus t = ±1. 

This implies that a^ 1 a 2 = p is a complex number, which has rational components. 
Then p 2 = (sis^ 1 ) 2 = z\z 2 x . Therefore z\z 2 = (pz 2 ) 2 , and Lemma [2.11 implies that 
pz 2 G G. Since Z\ and z 2 are squarefree, each Gaussian prime divisor of Z\ has 
multiplicity 1 in both z\ and z 2 , and the same holds for z 2 . Therefore z\ and z 2 are 
associates, and so p 2 is a unit in G. Thus N(p) = 1, and Lemma [2.11 shows that p is 
in G, so (3) of Theorem 15.41 is proved. 

Finally we prove (4). Suppose that (a, z) parameterizes the twin pair (6, rj). Then 
N(#) = N(a) 2 N(z) is a square if and only if N(z) is a square. Clearly if z G Z or 

G Z then N(z) is a square. 

Suppose that N(z) is a square and consider a Gaussian prime divisor 7r of z. As 2; 
is squarefree in G, the number 7r 2 does not divide z. Hence it = 1 + i is impossible, 
since all other Gaussian primes have odd norm, and so N(z) would be of the form 
4k + 2. Similarly, if p = N(n) is an odd prime (of the from 4k + 1), then p | N(z), 
and so the conjugate of 7r must also occur in z. Therefore z is indeed real or pure 
imaginary, and so the proof of Theorem 15.41 is complete. □ 
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Theorem 5.10. For a positive integer M denote by T(M) the number of twin pairs 
(9,rj) such that N(0) = N(r/) = M , and let a(s) be the sum of positive integer divisors 
of any integer s. Suppose that 

M = rtf...ifcqP ...gf , 

where Pi, ■ ■ ■ ,p m ar e primes = 1 (4) and q±, . . . , qe are primes = — 1 (4); we assume 
that all X r and [i s are positive. Then 

m i 

T(M) = 24n^)IIW S )' 

r=l s=l 

where 

g(p 2X ) = <V) + °(p x - 1 ) , g(p 2X+1 ) = Mp x ) , 

and 

h(q 2 ») = atf) + a^- 1 ) , h(q w ) = . 

In particular, T(M)/24 is a multiplicative function. If there exists a twin pair with 
norm M , then M is the sum of two squares. 

Proof. By Theorem 15.41 we have to count the number of pairs (a,z) GExG, where 
M = N(a) 2 N(z) and z is squarefree in G, and divide the number of solutions by 4 
because of (3). 

Writing; product of Gaussian primes we see that in the canonical form of 

N(z) every prime of the form 4k + 3 has exponent 2, every prime of the form 4k + 1 
has exponent 2 or 1, and the prime 2 has exponent 1 (or 0). If such a number 
t = N(z) is given, then the only freedom in determining z occurs at the primes p of 
the form 4k + 1. Indeed, p = niWi is a product of two Gaussian primes, and if the 
exponent of p in N(z) is 1, then we can decide whether to put 7r or W into z. We have 
to multiply the resulting z with the four Gaussian units. Thus if 4f(t) denotes the 
number of solutions for z with norm t, then / is a multiplicative function, which is 1 
or for every prime power except that f{p l ) = 2 when p = 1 (4). 

Corollary 12.141 allows us to count the number of integral quaternions a with given 
norm N(a). The result is 24 times a multiplicative function (the sum of odd divisors 
of N(a)). Now we go through all primes in the decomposition of M to see how we 
can split M into N(a) 2 N(z). 

If the prime is 2, then we must put 2 K into N(a) 2 if k is even, and must put 2 K_1 
into N(a) 2 if re is odd. By Corollary [2H we see that T(2 K ) = 24, and T(M) does 
not depend on re. 

Next consider a prime q r = — 1 (4). Then [i r must be even for a solution to exist, 
and we can either put the entire q^ r into N(a) 2 , or put q^ r ~ 2 into N(a) 2 and q r into z. 
This proves the formula in the theorem for h. 

Finally for p r = 1 (4) there are two cases to consider. If A r is even, then we can put 
or 2 copies of p r into N(z). If A r is odd, then we must put 1 copy of p r into N(z) 
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(and the corresponding Gaussian primes in z can be chosen in 2 ways). This proves 
the formula for g. 

Since we can put together the solutions for M from the solutions for the prime 
divisors of M independently, we get the formula in the theorem. □ 

Corollary 5.11. Let (v,w) be a pair of twins in Z 3 whose length is an integer. Then 
there exists an a £ E and an integer d such that the last two columns of dE(a) are 
either (v,w) or (—w,v). Therefore {v,w) can be extended to an icube. 

Proof. Let (y(v),V(w)) = (9,rj), parameterized by (a, z) where z is squarefree. By 
(4) of Theorem I5.4[ z is either real or pure imaginary, so z = d or z = di for some 
integer d. In the first case 9 = daja and 7] = dak a, so the last two columns of 
dE(a) are v and w. In the second case 9 = dak a and r\ = —daja, so the last two 
columns of dE(a) are — w and v. □ 

Corollary 5.12. // (u, v, w) is an icube, then there is an a EE, and d £ Z such that 
(u, v, w) and dE(a) can be obtained from each other by permuting and changing signs 
of columns. 

Proof. By Proposition 11.11 the edge length is an integer. Let a and d be given by 
Corollary 15.111 Then the columns of dE(a) and (±m, ±v, ±w) share two orthogonal 
vectors, so they share the third column of E(a) as well. □ 

Lemma 5.13. Suppose that a £ E, z is a squarefree Gaussian integer and 9 = azja 
is primitive. Then N(^) is the squarefree part of N(0). 

Proof. As N(0) = N(a) 2 N(z) it is sufficient to prove that N(z) is squarefree. Suppose 
that p 2 | N(z) for a prime < p £ Z. If p = 3 (4), then p is a Gaussian prime, so p \ z, 
contradicting the fact that 9 is primitive. If p = 2, then 2 | z, again a contradiction. 
Finally if p = ttW for some Gaussian prime tt, then 7r and W cannot both divide z, 
because then p would divide the primitive 9. On the other hand, the exponent of 
7r and 7f is at most 1 in z, since z is squarefree. Therefore N(z) cannot be divisible 
by p 2 , a contradiction. □ 

Proof of Theorem \1.5l We have proved the existence statement in Remark 14.31 Sup- 
pose that x is primitive. If (u, v, w) is an icube with edge length m such that 
x = au + bv + cw, then (u, v, w) is primitive, too. Therefore by Corollary 15.121 
(or by Corollary 13.91) we may assume that (u, v, w) = E(a) for some a £ E. Thus it 
is sufficient to deal with "Eulerian" cubic lattices. 

Suppose that x is contained in two such sublattices: V(x) = ai/3% 577 = a 2 /?2 c*2- By 
the uniqueness part of Theorem |4]2] there exists a unit e £ E such that 0J2 = a\e~ l and 
@2 and Pi are group-conjugates via e. Propositions 13. II and 13.21 show that the matrices 
E(ai) and E(a2) may differ only by permutations and sign changes of columns. 
Therefore the two cubic lattices are actually the same, proving the uniqueness part 
of the theorem. □ 
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Proof of Corollary \l.b\ We use the notation of the previous proof, let L denote the 
unique sublattice obtained there. Proposition 13.21 says that (3i and f3 2 have the same 
number of zero components. Therefore when considering the three cases of Corol- 
lary 11.61 (which are distinguished by the number of zero components of x relative 
to L), it does not matter which generating a we choose for L. 

We show that every twin of x is contained in L. Let 771 be a twin of 9 = V(x), and 
let (ai,zi) parameterize the pair (6,f]i) such that Z\ is squarefree. By Lemma [5.131 
N(zi) = n, hence N(«i) = m. Thus the sublattice generated by «i is L. Since Z\j 
has a zero component, in case (1) of Corollary 11.61 the vector x cannot have a twin. 

If the norm of x is a square, then 1 = n = N(zi), and therefore Z\ G {±1, ±2}. Thus 
x and each of its twins has exactly one nonzero component relative to L. Therefore x 
has exactly 4 twins. Conversely, if x has only one nonzero component relative to L, 
then its length is obviously an integer, since the same holds for the generating vectors 
of L. Hence (3) is proved. 

Finally suppose that none of the components of z\ is zero (so x has exactly one 
zero component). Let rj 2 be another twin of 9, parameterized by (a 2 ,z 2 ). Again, 
a 2 = a\e~ l and ezije' 1 = z 2 j for some unit e. However, not every unit e yields 
a twin of 9. Indeed, if e ^ Q = {±1, ±z, ±j, ±/c}, then Proposition 13.21 shows that 
group-conjugation by e induces a fixed point free permutation on the components of 
the vectors (while possibly changing some signs). We know that the first component 
of z\j and of z 2 j is zero. Therefore e £ Q can happen only if z%j and z 2 j have two 
nonzero components, which we have excluded. So e G Q. The two twins of 9 in 
question are rji = a\Z\k~ai and 

772 = a^klxi = —aiz^je^ie)^. 

Let us calculate this now. 

If e G {±1, ±«} then rj 2 = r}\. (This has been noted in (3) of Theorem 15.41 ) 
If e G {±j, ±fc}, then r] 2 = —rji (This is always obviously another twin of 9.) 
Thus if x has two nonzero components relative to L, then it has no more than two 

twins, completing the proof of Corollary 11.61 □ 

6. Twin-complete numbers 

In this section we first prove the characterization of twin-complete numbers given 
in Theorem 11.81 and then discuss Conjecture 11.91 

Lemma 6.1. If 4n is twin- complete, then so is n. 

Proof. Every pure quaternion whose norm is divisible by 4 is divisible by 2. This 
follows by looking at the coefficients mod 4 (or from Theorem 14. 6p . Thus if N(#) = n, 
then 29 has a twin 77, and so 77/2 is a twin of 9. □ 

Lemma 6.2. Let (3 G E be a primitive pure quaternion and m > an odd positive 
integer. Then there exists an a G E with norm m such that af3a is primitive. 
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Proof. It is sufficient to prove this when m is a prime, since we can go through the 
prime divisors of m one by one. Clearly (or by Theorem 14.61) . a/3 a is primitive if and 
only if it is not divisible by p = m. By Proposition 14. 7\ there is such an a (we have 
actually counted them). □ 

Proof of Theorem \1.8l Let n,m > be integers such that n is squarefree. First 
suppose that n is twin-complete. It is sufficient to prove that every primitive vector 
5 with norm nm 2 has a twin (since we can do induction on m) . By Theorem 14. 2\ 
5 = a/3 a for some a, (3 G E such that N(a) = m and N(/3) = n. Since n is twin- 
complete, (3 has a twin 7. We show that aja is a twin of 5, using Proposition 15.11 
Indeed, 

5a^a = a(3(aa)ja = maPja . 

This is a pure quaternion, since f3j is a pure quaternion. So we have proved one 
direction of the theorem. 

For the converse suppose that n > is square free and nm 2 is twin-complete. Then 
every vector (3 with norm n is primitive. By Lemma 16.11 we may assume that m is 
odd. For any given (3 Lemma [6.21 yields an a with norm m such that 9 = a/3 a is 
primitive. Since nm 2 is twin-complete, 9 has a twin. By Theorem I5.4[ this pair of 
twins can be parameterized by some (ckj, z) G E x G such that z is squarefree. Thus 
9 = aizjai, so by Lemma \5. 131 N(z) is the squarefree part of N(#), that is, N(z) = n 
and N(oti) = m. By the uniqueness statement of Theorem 14. 21 we get that (3 = ezje~ l 
for some unit e G E. But then ezke~ x is a twin of (3. Hence n is twin-complete. 

To prove the last statement of the theorem let n be square free. Clearly (a, b, 0) 
and (—6, a, 0) are twins, so if n is not the sum of three positive squares, but is the 
sum of at most two positive squares, then it is twin-complete. Conversely, suppose 
that n is twin-complete. Let (3 be a pure quaternion of norm n, we have to show 
that at least one of the three coordinates of the corresponding vector is zero. As 
n is twin-complete, (3 has a twin, so Theorem 15.41 implies that (3 = azj a for some 
(a,z) G E x G. Here n = N(/3) = N(a) 2 N(z), so N(a) = 1 and a is a unit. From 
Proposition 13.21 we get that at least one component of (3 is zero (since this is the case 
with zj). □ 

Now we discuss Conjecture 11.91 Let 

S D {1, 2, 5, 10, 13, 37, 58, 85, 130} 

denote the list of those squarefree numbers that can be written as a sum of at most 
two positive squares, but not as a sum of three positive squares. It has been known 
since |CCC59j that this list is finite, and if the conjecture fails, there is at most one 
number in S not listed above ( [ Wei73] . |Gro85j ). 

In |Mor60] it is shown that for an integer n G S, the only nonnegative solutions of 



xy + yz + zx = n 
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when n = 2 (4) are given by xyz = 0, and when n = 1 (4) are given by either xyz = 
or x = d, y = d, z = (n — d 2 )/2d, where d is any divisor of n with d 2 < n. Either way, 
for such numbers n the above equation has no solution with three distinct positive 
integers x, y, z. 

This characterization of the numbers in 5* allows us to see the relationship they 
bear with Euler's numeri idonei. Euler defined a numerus idoneus to be an integer N 
such that, for any positive integer m, if 

m = x 2 ±Ny 2 , (x 2 ,Ny 2 ) = l, x,y>0 

has a unique solution, then m is of the form 2 a p k , a G {0, 1}, k > 1, p is a prime. 

Euler was aware of 65 numeri idonei, and it is widely believed and conjectured 
that his list is complete ( |Rib00] ) . It was S. Chowla who proved that there are only 
finitely many numeri idonei ( [Cho34j ). P. J. Weinberger improved Chowla's result to 
show that if Euler's list is not complete, then at most one numerus idoneus is absent 
and it must be greater than 2 • 10 11 ( |Wei73j ). 

Using Theorem 3.22 of [Cox89j it can be shown that an integer N is a numerus 
idoneus if and only if it cannot be expressed as xy + yz + zx with < x < y < z. 
Combining this with the characterization by |Mor60j described above, we see that 
every integer in S is also one of Euler's numeri idonei. Checking Euler's list of the 
65 numeri idonei (the greatest of which is only 1848) against the properties listed in 
Conjecture 11.91 one sees that, indeed, Conjecture II .91 is true if Euler's list is complete. 

If we only consider those integers n for which there is no representation of the 
form xy + yz + zx with 1 < x < y < z, i.e. the even, squarefree, twin-complete 
numbers, then we have from [BCOOj that such an n can only be absent from the list 
of Conjecture 11.91 if the Generalized Riemann Hypothesis fails. 

7. Geometric considerations 

Our goal in this section is to shed light on the main results of the paper by using 
elementary methods. Doing so, we reprove several results in a geometrical way. The 
vectors that we work with will be assumed to be non-zero, unless otherwise stated. 

Our first result strengthens Lemma 15.51 in the case of equally long vectors. 

Proposition 7.1. Let x,y E Z 3 be two orthogonal vectors of equal norm (length 
squared) N(x) = N(y) = nm 2 , where n,m G N ; n is squarefree. We claim that the 
cross product x x y is divisible in Z 3 by nm. 

Proof. Let x = (xi,X2,xz) and y = (2/1, 2/2, 2/3) - Our hypotheses mean that 

x\ + x\ + x\ = nm 2 = yf + yj + y\, 
XiVi + x 2 V2 + x 3 y 3 = 0. 
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From the second equation we get that x\y\ = x\y\ + x\y\ + 2x\x 2 y\y 2 . We compute 
the square of the third coordinate X\y 2 — x 2 y\ of x x y as follows: 

(x x y 2 - x 2 yif = (x\ + xl) (y\ + yl) - x\y\ 

= (nm 2 - x 2 3 )(nm 2 - yl) - x\y\ 

2 / 2 2 2\ 

= nm (nm —x 3 — y s ). 

As n is squarefree, n 2 m 2 divides the square (xiy 2 — x 2 yi) 2 , hence nm divides the third 
coordinate of x x y. By a similar argument it divides all three coordinates. □ 

Corollary 7.2 (The case n = 1, reproving 15. lip . If x,y G Z 3 have length m e Z and 
ili/, i/ien x x y is divisible (in Z 3 j fry m, hence the vector (x x y)/m extends the 
twin pair (x, y) to an icube. □ 

Corollary 7.3 (The case m = 1). Ifx,ye Z 3 where N(x) = N(y) = n is squarefree 
and x J_ y, then x xy is divisible by n. Hence (x xy)/n is one of the six unit vectors 
ofl?. In particular, x and y lie in one of the three coordinate planes. □ 

Next we reprove the case of Theorem II .51 for vectors of integer length by elementary 
means. This also implies Theorem [LH By a "square sublattice" we mean a sublattice 
generated by a twin pair. Sublattices are subgroups of Z 3 , and therefore it makes 
sense to speak about the index of them in the sense of group theory. 

Theorem 7.4. Let v = (a, b, c) be a primitive vector in Z 3 with integral length £, 
and let L = {x £ Z 3 | x _L v}. Then amongst the square sublattices of L there is a 
greatest one K (with respect to set inclusion), and the index i of K in L equals £. 

Proof. A non-trivial sublattice of L is square if and only if it is invariant under the 
90° rotation R about the vector v. Since the sublattice generated by all i?-invariant 
sublattices of L is again /^-invariant, there is a greatest such lattice K. 

The vectors r = (0, —cl, bl) and s = (b 2 + c 2 , — ab, —ac) in L are of equal length 
and orthogonal to each other, thus they span a square sublattice of L. The area of 
the square spanned by the vectors r and s is £ 2 (b 2 + c 2 ). Denote by K\ the lattice 
generated by these vectors. By analogous constructions we obtain two other square 
sublattices K 2 and K%, in which the area of the fundamental square is £ 2 (c 2 + a 2 ) and 
I 2 (a 2 + b 2 ), respectively. The area of a fundamental parallelogram of L is well-known 
to be £, thus the indices of the lattices Ki in L are £(b 2 + c 2 ), £(c 2 + a 2 ), and £(a 2 + b 2 ). 
Our hypothesis gcd(a, 6, c) = 1 easily implies that gcd(6 2 + c 2 , c 2 + a 2 , a 2 + b 2 ) = 1. 
(Here we take advantage of the fact that out of the numbers a, b, and c exactly one 
is odd, since the sum of their squares is a square.) This, in turn, implies that the 
index i of K in L divides the number £ = gcd (£(b 2 + c 2 ),£(c 2 + a 2 ),£(a 2 + b 2 )) . 

Let the vectors u and w span a fundamental square (of area i£) of the lattice K. 
According to Proposition I7.1[ the length-square i£ of the vectors u and w divides the 
square of each coordinate of the vector u x w = iv, that is, £ \ ia 2 , £ \ ib 2 , and £ \ ic 2 , 
hence £ \ gcd(ia 2 , ib 2 , ic 2 ) = i, proving the claim i = £. □ 
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Corollary 7.5. Using the hypotheses of the previous theorem, there exists exactly 
one cubic sublattice C of Z 3 containing v as a shortest nonzero vector. 

Proof. The greatest square sublattice K of L has index I in L, hence its shortest 
vectors share the length I with v, consequently K and the vector v span the required 
cubic sublattice. 

On the other hand, if C is a cubic lattice containing v as a shortest (non-zero) 
vector, then the orthocomplement lattice K' of v in C is a square sublattice of L 
with index £, hence K' coincides with the greatest square sublattice K of L. □ 

A corollary of Theorem 15.101 is that the norm of twins is the sum of two squares. 
Here we prove this assertion by elementary means. 

Proposition 7.6. The norm of twin vectors in Z 3 is the sum of two squares. 

Proof. Let x, y E Z 3 be twin vectors whose norm is M. Consider the lattice L of 
all integer vectors z E Z 3 lying in the plane spanned by x and y, and let R be the 
rotation through 90° in the plane of L. Set V = LC)R(L). Plainly, L' is a non-trivial 
lattice invariant under the rotation R. Therefore, a shortest non-zero element z E V 
and Rz together span V . Let I be the length of z. The vector x x y has integer 
length M, hence it has a twin w in V by Corollary 17.51 Since x,w E V and L' is 
spanned by z and Rz, there are integers a,b,c,d such that \fM = £y/ a 2 + b 2 and 
M = £y/c 2 + d 2 . Dividing these equations and taking squares we get that the integer 
M is equal to (c 2 + d 2 ) / (a 2 + b 2 ), hence, each of its prime factors of form 4k + 3 occur 
on an even power. □ 



This section claims explicitly the characterization of Pythagorean quadruples, lists 
the theorems of the paper from which this result follows, and gives another, short, 
constructive proof using Gaussian integers. Consider the Diophantine equation 



It is sufficient to find the primitive solutions, that is, when gcd(a, b,c,d) = 1. Consid- 
ering this quadratic equation mod 4, we realize that d is odd, and among the integers 
a, b, c exactly one is odd. We may assume that this number is a. 

Theorem 8.1. The quadruple (a,b,c,d), where a 2 + b 2 + c 2 = d 2 , gcd(a, b, c) = 1, 
a is odd, d > 0, possesses exactly four Euler parameterizations 



8. Pythagorean quadruples 



2 i 1.2 , 2 j2 

a + b + c = d . 



d 



a 



b 



c 
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with some integers m, n, p, q. By setting n = q = we get the usual parameterization 
of Pythagorean triples (where b = 0). 

This result is found as formula (11) on page 31 of [Carl5j . It obviously follows from 
the fact that (a,b,c) is contained in a unique cubic lattice (see Theorem II .51) . which 
has been proved using quaternions in Section [5l and also by geometrical methods in 
Corollary 17.51 

Lemma 8.2. Let x,y G Z be positive integers. Suppose that xy = 77 with 7 G G 
and gcd(x, y, 7, 7) = 1. Then there exist integers m, n,p, q G Z such that 

x = (m + ni) (m — ni), 7 = (m + ni) (p + qi) , 

V — (p + Q.i){p — qi), 7 = (m — ni){p — qi) . 

Every such quadruple (m,n,p,q) automatically satisfies that 

m + ni = gcd(x, 7), 

p + qi = gcd(y,7). 

Furthermore, the exact number of such quadruples is four. 

Proof. We use the notation gcd(a,/3) for Gaussian integers, too. Of course this 
number is unique only up to multiplication by a unit element of the ring G. We use 
~ to denote associates. 

Let 6 = gcd(x, 7) = m + ni and r] = gcd(y, 7) = p + qi. Taking conjugates we get 
gcd(x, 7) = 9 and gcd(y,7~) = rj. Clearly gcd(#, rj) \ gcd(x, y, 7, 7) = 1. From 

9 rj 9 rj 

we see that x/9 and /rj divide each other, since gcd(x/6 l , 7/6 1 ) = gcd(y/fj,'y/fj) = 1. 

Now gcd(6 ) , rj) = 1, hence j/rj ~ gcd(9'j/fj,rj' : y/fj) ~ gcd(x, 7) = 9. Modifying 9 
by multiplying it with a suitable Gaussian unit we get 7/7/ = 9. Then 7 = 9rj and 
x ~ #7/77 = 99. Both are positive integers, so x = 99. Hence y = 77/x = rffj. 

Since 9 = m + ni is uniquely determined after fixing one of the four possible 
values of rj = gcd(y, 7), we get that the number of all quadruples (m, n,p, q) with the 
required properties is exactly four, thus completing the proof of the lemma. □ 

Proof of Theorem \8.1[ Write the Pythagorean quadruple equation as 

c\ 2 ( b\ 2 d + a d — a 



< 8 - 3 > UJ + V2 
or, equivalent ly, 

. „. f c b\ ( c b \ d + a d — a 

in G. 



2 2 / V 2 2 
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Apply Lemma for x = (d + a)/2, y = (d - a)/2 and 7 = -c/2 + (b/2)i. We 
must show that gcd((d + a) /2, (d - a) /2, -c/2 + (6/2)i, -c/2 - (b/2)i) = 1. 

Suppose that £ G G divides (d — a)/2 and (d + a)/2. It must then divide a and d. 
Suppose further that 5 also divides —c/2 + (b/2)i and —c/2 — (b/2)i. It must then 
divide their sum, which is — c, and their difference, which is hi. Therefore, 5 divides 
a, b and c. By assumption, gcd(a, b, c) = 1, and so we must have 5 G {±1, ±z}. 

We may now apply Lemma [8.21 to yield integers m, n,p, q such that 

= (m + ni)(m — ni), 



2 

d — a 



(p + qi){p- qi), 



~ + ^* = (m + ni)(p + qi). 

The above equations are easily seen to be equivalent to the conclusion of the theorem. 
The number of such Euler parameterizations (m,n,p,q) is exactly four according to 
the last statement of the previous lemma. □ 

Remark 8.5. We see from Lemma 18.21 and Theorem 18.11 that in order to determine 
a set of parameters m, n,p, q for a given quadruple (a, b, c, d), we need only compute 
gcd ((d + a)/2, — c/2 + (b/2)i) to yield m + ni, and then divide — c/2 + (b/2)i by this 
result to obtain p + qi. If b = 0, then m + ni and p + qi can be chosen to be real. 

9. Problems and conjectures 

Problem 9.1 (Main problem). Investigate construction, counting, extension for 
higher dimensional icubes. Describe the possible edge lengths of icubes of a given 
dimension k in Z n , and twin-completeness. Generalize the results in this paper. 

We make some remarks and state two conjectures. First note that if the dimension 
is even, then every vector has a twin: 

(a, b, ... , c, d) and (—6, a, . . . , —d, c) 

are twins. Therefore odd dimensions seem more interesting. 

Call an integral vector odd if every component is odd. Looking at the vectors 
mod 2 we see that if the dimension is odd, then odd vectors do not have a twin. 

Conjecture 9.2. Suppose that the dimension n > 5 is odd. Then every non-odd 
vector has a twin. 

Every odd vector of length y/d satisfies that d = n (8). Therefore the following 
conjecture is weaker. 

Conjecture 9.3. Suppose that the dimension n > 5 is odd and d ^ n (8). Then 
every vector of length \fd has a twin. 
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Since the number of vectors of a given length grows rapidly with the dimension, it 
may be possible to use analytic methods in investigating these questions. 
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